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ABSTRACT 


A  review  of  the  current  theory  and  methods  for  the  computa¬ 
tion  of  maximum  flow  In  networks  Is  presented  along  with  a  simplified 
method  for  determination  of  a  feasible  -Mow  In  networks  with  upper 
and  lower  bounded  arcs.  A  computational  procedure  Is  presented  which 
Is  used  to  calculate  the  maximum  flow  for  a  general  network.  Pie 
network  Is  reduced  to  an  equivalent  basic  network.  An  associated 
network  Is  used  to  compute  a  feasible,  then  the  maximum  flow  for 
the  basic  network.  A  computer  program  Is  Included  for  use  In 
computation  of  maximal  flows  In  large  networks. 
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CHAPTER  I 


INTRODUCTION 

The  growing  uses  for  mathematical  programming  models  In  trans 
portation  and  communication  networks  have  prompted  study  in  the  devel 
opment  of  computational  methods  to  determine  maximal  network  flows. 
The  subject  of  this  thesis  is  the  mathematical  problem  of  determining 
maximal  steady  state  flows  in  networks  which  are  subject  to  capacity 
limitations  on  the  arcs  and  nodes  of  the  network. 

An  efficient  and  widely  used  method  to  compute  maximal  flow 
was  developed  by  Ford  and  Fulkerson  [7].  This  method  along  with  a 
review  of  the  current  theory  will  be  presented  In  Chapter  II  In 
order  to  give  thevreader  the  theoretical  foundation  upon  which  the 
remainder  of  the  thesis  is  based.  Chapter  III  Introduces  a  new 
method  which  utilizes  the  Ford  and  Fulkerson  algorithm  to  compute  the 
maximal  flow  In  networks  with  non-zero  lower  bound  limitations  on  the 
flow  In  the  arcs.  Since  many  networks  of  Interest  are  more  general 
In  nature,  the  geierallzed  network  Is  introduced  In  Chapter  IV  along 
with  a  computational  procedure  for  the  determination  of  maximal  flow. 
To  Illustrate  the  techniques  proposed,  an  example  Is  Included.  The 
procedure  can  be  used  to  solve  smal  l  problems  by  hand,  but  the  amount 
of  work  Increases  rapidly  with  the  size  of  the  network.  Therefore,  a 
computer  program  which  may  be  used  to  solve  larger  problems  Is 
included  In  the  appendices.  A  brief  description  of  the  program  Is 
found  In  Chapter  IV. 
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Definitions  and  Symbols 

In  order  to  establish  a  common  understanding  as  to  the  exact 
meaning  of  various  terms  and  symbols  used  to  present  the  material  of 
this  thesis,  an  Initial  set  of  definitions  Is  presented.  Other  terms 
will  be  defined  throughout  the  material  as  needed  to  facilitate  the 
presentation. 

Define  a  network ,  G,  as  a  graph  which  can  be  represented  on 
a  plane  In  such  a  way  that  the  set,  N,  of  nodes  are  distinct  points 
and  the  set,  A,  of  arcs  are  simple  curves  which  connect  two  distinct 
nodes.  Furthermore,  no  two  arcs  can  meet  except  at  the  nodes  which 
are  their  extremities.  The  graph  contains  no  loops. 

Define  (o~(N)  to  be  the  set  of  arcs  or  flows  Incident  to  a  set 
of  nodes,  N,  and  (/(N)  to  be  the  set  of  arcs  or  flows  Incident  from  N. 

A  directed  arc,  ai  y  Is  defined  as  an  arc  Incident  from  a  node 
1  and  Incident  to  a  node  J. 

An  undirected  arc  conr.?rts  node  1  and  node  j  without  having  an 
orientation. 

Flow  can  be  thought  of  as  a  value  of  the  steady  state  rate  of 
movement  of  a  homogeneous  commodity  along  a  path  or  channel.  Positive 
flow  In  a  directed  arc,  a^  j,  will  move  with  the  orientation  of  the 
arc  and  will  be  denoted  by  the  symbol,  y^  j.  Positive  flow  In  an 
undirected  arc,  denoted  y^  y  may  move  In  either  direction  but  not 
In  opposite  directions  at  the  same  time. 

To  each  arc,  a^  j  c  A,  there  will  be  associated  two  real 
numbers,  b^  j,  and  c^  j,  which  represent,  respectfully,  the  minimum 
and  maximum  allowable  flow  In  that  arc.  The  bounds  associated  with 
undirected  arcs,  a^  j,  are  denoted  b^  j,  and  c^  y 
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An  arc  Is  saturated  If  y^  j  ■  j. 

Define  Y.  as  the  set  of  flows  in  a  network  G. 

Initially,  the  capacities  of  the  nodes  will  be  assumed  to  be 
Infinite.  The  nodes  will  be  classified  according  to  the  value  of  the 
difference  of  the  sum  of  the  flows  Into  the  node  minus  the  value  of 
the  sum  of  the  flows  out  of  the  node.  A  node  Is  designated  as  a 
source  node  If  this  difference  Is  negative,  and  as  a  sink  node  If  this 
difference  Is  positive.  For  all  Intermediary  nodes,  the  sum  of  the 
flows  Into  a  node  must  equal  the  sum  of  the  flows  out.  The  set  of 
source.  Intermediary  and  sink  nodes  will  be  denoted  as  S,  I,  and  T 
respectfully. 

Define  F  as  the  value  of  the  flow  In  a  network. 

F  -  i  y{  t 

yi,t» <T)  ■ 

A  flow  Is  feasible  If  and  only  if: 


(1)  b1,J -y1,J -C1,J  Va1,j£A 

(2)  F  -  z  +  ys  , 

yStjcu> \s)  S»J 


1  -  t  '  F 

yuto»  (T)  1*t 

Z  y  -  E 

y1  >jCtD+CJ) 


Define  a  chain  as  a  sequence  of  arcs  such  that  each  arc,  a.  . 

■  iJ 

1,J  c  I,  Is  connected  to  an  arc  a^  1  at  node  1  and  to  arc  a^  1  at  node  j. 
A  cycle  Is  a  sequence  of  arcs  where: 

(1)  The  sequence  Is  a  chain. 

(2)  The  sequence  does  not  use  the  same  arc  twice. 

(3)  The  Initial  and  terminal  nodes  of  the  chain  coincide. 


A  cut  Is  defined  as  any  set  of  directed  arcs  containing  at 
least  one  arc  from  every  chain  of  positive  capacity  joining  the 
source  to  the  sink.  The  value  of  the  cut  Is  equal  to  the  sum  of  the 
capacities  of  the  arcs  of  the  cuts. 

Basic  Network 

Much  of  the  Initial  development  of  the  material  In  this  thesis 
deals  with  computation  of  flows  In  simple  or  basic  networks.  There- 
fore,  for  the  purpose  of  this  paper,  a  basic  network,  G,  Is  defined 
as  a  network  which  has  the  following  properties. 

1.  The  network  has  a  return  arc,  aR  y  which  Is  added  to  the 
network  for  computational  purposes  only. 

2.  Node  1  has  only  one  arc,  anj,  Incident  to  it.  Node  1  Is 
designated  as  the  source.  In  the  basic  network  there  Is 
only  one  source  node. 

3.  Node  n  has  only  one  arc,  an^,  Incident  from  It.  Node  n 
Is  designated  as  the  sink.  In  the  basic  network  there  Is 
only  one  sink  node, 

4.  For  each  arc,  a^  j,  tnc  associated  c^  j ,  Is  positive  and 

b1,J  "  °* 

5.  Two  distinct  nodes  (1,j)  may  be  connected  by  only  one 
directed  arc.  The  network  contains  no  undirected  arcs. 


CHAPTER  II 


REVIEW  OF  THE  LITERATURE 

In  order  to  acquaint  the  reader  with  the  subject  of  maximal 
flow  In  networks,  a  review  of  the  current  theory  and  primary  coafclna 
tlonal  techniques  is  presented  In  this  chapter. 


Basic  Network  Theory 

A  great  deal  of  theoretical  work  In  networks  has  already  been 
dore.  A  few  of  the  relevant  theorems  are  presented  here  In  order  to 


provide  justification  and  a  rationale  for  the  combinatorial  techniques 
and  algo.lthms  that  are  presented  In  later  sections. 

Theorem  1.  Opposite  directed  flows  on  the  same  arc  cancel  [4]. 
This  theorem  means  that  given  an  undirected  arc,  a^  y  with  flows 
y^  j  >  0  and  y^  ^  >.  0,  these  two  flows  can  be  replaced  by  the  flows 

y'i.j  and  y‘j.i 


Theorem  2.  A  set  Y  with  F>0,  satisfying  the  capacity  con¬ 
straints  and  node  conservation  equations  can  be  decomposed  Into  a  set 
of  positive  chain  flows  from  the  source  to  the  sink  and  a  set  of 
circular  flows  such  that  the  direction  of  positive  flows  In  any 
common  arc  Is  the  same  for  all  chains  [4]. 

Theorem  3.  There  exists  a  positive  flow  from  the  source  to 
the  sink  If  there  exists  at  least  one  chain  of  arcs  with  positive 
capacity  from  the  source  to  the  sink  [4], 

Theorem  4.  A  flow  F°  Is  maximal  If  and  only  If  the  maximal 
flow  Is  zero  In  a  second  network  formed  by  replacing  c4  4  by  c4  ,  - 

■  »  J  *  •  J 

y1.J  V  a1,J- A* 
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Theorem  5.  Given  any  partition  of  nodes  into  two  sets,  where 
the  first  set  Includes  the  source  node  and  the  second  set  the  sink, 
then  a  feasible  solution  F  Is  maximal,  if  every  arc,  a^  j,  that  Joins 
a  node  In  the  first  set  to  a  node  of  the  second  Is  saturated  [4]. 

Theorem  6.  The  maximum  flow  In  a  network  Is  equal  to  the 
minimum  cut. 

Tree  Method 

Dantzlg  [4]  has  developed  a  systematic  procedure  for  com¬ 
puting  the  maximal  flow  In  a  network  by  using  chains  connecting  the 
source  and  the  sink  of  the  network.  A  knowledge  of  the  tree  method 
will  give  a  better  understanding  of  the  Ford  and  Fulkerson  algorithm 
which  Is  presented  In  the  next  section  because  the  two  methods  are 
Identical  In  principle. 

A  tree  [1]  Is  defined  as  a  connected  graph  with  at  least  two 
nodes  and  no  cycles. 

Consider  a  basic  network  such  that: 

bi.J  *  0  ’  C1  .J  -  0  V  a1.J  -A 

The  arcs  may  be  directed  or  undirected.  Initially  all  arcs 
are  unmarked  and  the  feasible  flow  Is  zero.  In  order  to  determine 
the  maximum  flow: 

1.  Create  a  tree  such  that  the  arcs  are  unsaturated. 

2.  Select  two  sub  trees  -  one  Tj,  branching  out  from  the  source, 
node  1,  and  the  other,  Tn,  branching  out  from  the  sink,  n,  such  that 
each  node  Is  reached  by  just  one  arc  of  the  tree. 

3.  Choose  one  unsaturated  arc  which  connects  the  two  trees.  Thus, 
there  will  be  just  one  chain  from  node  1  to  node  n.  If  none  can  be 
found,  go  to  step  7. 
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4.  Determine  the  amount  that  the  flow  can  be  Increased  along  the 
chain  from  node  1  to  node  n,  by  the  following  procedure.  Determine 
the  amount  that  the  flow  can  be  Increased,  AX,  where  AX  -  min  (c^  j- 
y1  j*  l)’  Flow  a^on9  the  chain  directed  toward  the  sink  Is 
designated  j  and  flow  directed  toward  the  source  along  the  chain 
is  designated  y^  y  There  are  two  cases: 

(a)  AX  >  0,  Increase  the  flow  by  AX  along  the  chain,  go  to 
next  step. 

(b)  AX  =  0.  This  means  one  arc  In  the  chain  Is  saturated 
by  a  flow  directed  toward  the  sink,  go  to  next  step. 

5.  Mark  one  of  the  saturated  arcs. 

6.  Eliminate  the  saturated  arc  from  the  chain.  Thus,  there  are  two 
sub- trees.  Return  to  step  3. 

7.  Check  all  marked  arcs  joining  Tj  to  T  .  If  all  these  arcs  have 
flow  from  T-j  to  Tn  then  the  optimum  has  been  reached.  If  any  arc 
has  flow  from  Tn  to  T-j ,  use  this  arc  to  connect  T^  to  Tn  and  proceed 
with  step  4. 

Ford  and  Fulkerson  Algorithm 

Given  the  basic  network,  the  most  direct  method  developed  to 
date  of  computing  the  maximum  flow  Is  the  algorithm  developed  by  Ford 
and  Fulkerson.  This  method  Is  the  one  most  cited  In  the  literature 
concerning  maximal  network  flow.  The  procedure  Is  outlined  below. 

1.  Begin  with  a  feasible  flow  In  0  In  all  arcs. 

2.  Label  each  node  as  follows: 

a.  Label  node  1,  with  the  number  1. 


■p 


i 

I  i 


b.  If  node  i  Is  labeled  and  node  j  Is  not  labeled,  label 

node  j  with  the  number  (1)  If: 

(1)  a.|  j  exists  and  y-|  j  <  j 

(2)  aj  .  exists  and  yj  ^  >  0 

c.  If  the  sink,  node  n,  becomes  labeled  by  this  procedure, 
then  the  flow  from  node  1  to  node  n  can  be  increased. 

3.  If  node  n  is  labeled,  construct  a  simple  chain  from  node  1 
to  node  n  by  back  tracking  from  n  to  1  using  the  labels 

on  the  nodes. 

a.  The  chain  created  will  have  the  property  that  all  arcs, 
a,j  j  on  the  chain  directed  toward  the  sink  have  a  flow 
0  £  y^  j  1  j  and  all  arcs  a-|  k  directed  toward  the 
source  will  have  a  flow  y-j  k  >  0. 

b.  The  flow  can  be  Increased  by  an  amount  AX  where: 

AX  =  m1n(c1  -  y1  j,  y1>k)  V  a^j,  a1>k  In  chain. 

4.  Increase  the  flow  along  the  chain  by  an  amount  AX  and 

return  to  step  2. 

5.  This  procedure  Is  repeated  until  the  sink  cannot  be 
labeled.  Flow  Is  maximized  when  this  Is  the  case. 

Now  consider  the  more  general  problem  In  which  the  lower  bounds 

on  the  directed  arcs  are  no  longer  assumed  to  be  zero.  Thus  the  flow 

In  each  arc  Is  bounded  both  above  and  below. 


Such  that:  -  ®  £  b.j  j£ci  j£+” 

Theorem  7r  In  the  case  of  networks  with  lower  bounds  on  the 


a1,j 


arcs,  C  is  defined  as  the  subset  of  N  such  that  if  node  1  ^  C  the  node 
n  ^  C  or  If  node  1  £  C  then  node  n  £  C. 


A  flow  y.  .  exists  and  is  feasible  If  and  only  If: 

'  »J 

b1,j  lyiJ  £C1J  V  a1.j  -  A 

2.  Z  c»  i  ^  Z  .  b.  , 

»1  j  c  o)-(C)  a1  j  c  a)  (C)  1,J 

This  more  general  type  of  problem  can  be  satisfactorily  worked 
by  using  the  Ford  and  Fulkerson  algorithm  if  a  feasible  flow  can  be 
found. 

The  Ford  and  Fulkerson  algorithm  can  be  used  to  solve  for  a 
maximum  flow  In  an  associated  network  G'  which  Is  derived  from  the 
original  network,  G.  This  method  Is  outlined  below. 

1.  Create  an  associated  network  G'  such  that  G'  contains 
the  nodes  and  arcs  of  the  original  network  G  with  arc 
capacities  c'i  j  a  ci  j  “  bi  j  •  The  bound, 
b'.j  j  equals  zero. 

2.  Add  a  new  source,  labeled  0,  and  a  new  sink,  labeled  m, 

(m  =  n  +  1),  to  G'. 

a.  If  the  arc,  a^  j,  exists,  then  construct  two  new 
arcs  according  to  the  value  of  the  lower  bound,  b.  .. 

If  bi  j  >  °»  then  construct  the  arcs,  a'^  m  of  capacity 

c'l.m  ■  b1.j  a"a  *'o.j  wUh  “P*014*  c'o.J  *  b1,j  • 

If  b.  .  <0,  then  construct  the  arcs,  a'  7>#1th 

» Of! 

capacity  c'  .  =  -  b.  .  and  a' .  with  capacity 

o  >  i  1 1 J  J  >ni 

C1 ,  m  =  4  ■ 

j  ,m  1 ,  J 

b.  The  Ford  and  Fulkerson  algorithm  can  then  be  used  to 
determine  the  maximum  flow  In  G'  from  node  0  to  node  m. 

3.  Two  cases  can  result  from  step  2b. 

a.  Case  one:  If  the  flow  In  each  arc,  a'  .  and  a'. 

o,J  J,m 
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j  c  N',  equals  the  capacity  of  each  arc,  then  the  flow 
Is  feasible  for  the  network  6  and  the  problem  has  a 
solution. 

b.  Case  two:  If  for  any  arc,  a'  <  or  a*.  „  j  c  A*, 

o,j  j»m  °  - 

the  flow  In  that  arc  does  not  equal  the  capacity  of 
the  arc,  the  problem  has  no  solution. 

4.  If  a  feasible  flow  for  G  can  be  found  by  the  above  pro¬ 
cedure,  transfer  the  flow  from  G'  to  G  by  using  the 
following  transformation: 

y1.J*y'lJ  +  bl.J  Vsi.j£A 

5.  Ford  and  Fulkerson's  algorithm  Is  then  used  to  determine 

the  maximum  flow  using  the  feasible  flow  as  a  starting 
basis. 
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CHAPTER  III 


A  SIMPLIFIED  METHOD  FOR  COMPUTATION  OF  MAXIMAL  FLOW 
To  compute  a  maximum  flow  In  networks  with  non-zero  lower 
bounded  arcs,  an  Initial  feasible  flow  must  be  determined.  The  method 
presented  In  the  last  chapter  produces  a  complex  associated  network, 
G',  aven  with  a  relatively  simple  network,  6.  The  principle  advantage 
of  the  method  presented  here  [10]  Is  that  an  initial  feasible  flow 
can  be  determined  by  use  of  a  more  simple  network,  G',  than  in  the 
previous  method. 


Proposed  Algorithm 

The  development  of  the  algorithm  Is  based  upon  a  consequence 
of  the  flow  balance  equations.  In  order  to  provide  a  rationale  for 
the  procedures  presented,  the  development  of  the  associated  network 
will  be  presented  now. 

Given  a  basic  network  G  where  the  flow  In  each  arc  Is  bounded 
as  follows: 

'  " - bi .j  iyi,j 1  ci,j +  *  *t.j  £  A 
The  problem  of  determining  the  maximum  flow  in  G  can  be  stated  as 

follows: 


Maximize: 
Subject  to: 


Vi 
■yj.1 


E  Vi  1  -  l  +  yH»0  VIcN 

yjj  £  w  (O  y1fj  £ w  (D 


Vai,j^A 


bi,j  lyu  £Ci.J 

Now  transform  the  entire  network  G  Into  an  associated  network 

G*  using  the  change  of  variables: 


y1,j  =  y'i,j  +  b1,j 


Va1fj  cA 
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ru.  'vw  d  /■ 


V^jWS 


the  problem  then  becomes: 
Maximize: 

Subject  to: 


1  y'i  i  •  1  +  i +  -  bn 

y'uc.'H!  J’1  y'^cuO)  1,J  c  u  (1)  J>1 

E  .  b,  .  ■  0  V  1  C  N 

b,,j  i> (1>  Ui 


and:  0i*'l>jic'l>j'c1>j-b1.j 


i-cw  u  u,.  £  i  b.  j  =  dj  i  u  ii 

bj  1  c  uf(1)  J>1  b1  j  c  w  (1)  1,J  1 

Note  two  things.  First  the  sum,  E  d. ,  Is  equal  to  zero. 

1c  N  1 

Second,  If  a  flow  rf  -d^  was  passed  through  each  node,  the  problem 
would  be  exactly  the  problem  of  flow  maximization  In  a  basic  network. 
Hence,  a  technique  similar  to  Berge's  [2]  can  be  used  In  which  new 
nodes  0  and  m  *  n+1  are  added  to  the  network  G'.  Therefore,  only  a 
maximum  of  one  arc  per  node  Is  added  to  G'  In  this  method. 


1  c  N 


zero. 


Procedure 

Using  the  results  of  the  previous  section  the  solution  to  the 
maximal  flow  In  the  above  network  can  be  determined  using  the  following 
procedure. 

1.  Add  a  new  source  and  sink  node,  denoted  0  and  m  (m  -  n  +  1), 
respectfully  to  G'. 

2.  The  arcs  a'  t  and  a'.,  „  are  added  as  follows.  If  d.  >  0, 

o,l  j,m  1  * 

theflan  arc,  a'Q  ^ ,  is  constructed  with  capacity  d^.  If  d^  <  0,  an 
arc,  a'j  m,  Is  constructed  with  capacity  -d^.  If  d^  e  0,  no  arc  Is 
added  to  G'  for  node  1. 


3.  The  Ford  and  Fulkerson  algorithm  is  then  used  to  maximize 
flow  In  6'  from  node  0  to  node  m.  Again  as  in  the  previous  method 
this  maximum  flow  is  equivalent  to  a  feasible  flow  in  the  network  G 
If  and  only  If: 

y*0,1  °  d1  ¥,3d>0 

V  i  3  d  <  0 

4.  If  the  problem  Is  solvable  continue  by  eliminating  all 
arcs  Incident  from  node  0  and  Incident  to  node  m.  Node  1  and  node  n 
are  now  the  source  and  sink  respectfully  for  G'. 

5.  Using  the  flow  generated  in  step  3  as  a  starting  basis, 
maximize  flow  from  node  n  to  node  1  by  again  using  the  Ford  and 
Fulkerson  Algorithm. 

6.  Upon  completion  of  step  5,  transfer  flow  from  G'  to  G  by 
use  of  the  transformation  equation: 

yi.j +  buj  Vai.j^A 

The  maximum  flow  In  G  Is  equal  to  yn  -j. 
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CHAPTER  IV 


PROCEDURE  FOR  THE  SOLUTION  TO  THE  GENERALIZED  MAXIMUM  FLOW  PROBLEM 

A  proposed  methology  for  the  computation  of  a  solution  to  the 
generalized  maximal  flow  type  problem  as  defined  below  Is  presented. 
The  generalized  network  is  first  reduced  to  a  basic  network  and  then 
the  technique  presented  In  Chapter  III  Is  used  to  solve  for  the  equlv*- 
talent  maximum  flow  In  the  basic  network.  An  example  Is  presented  in 
order  to  demonstrate  the  techniques  discussed  in  this  chapter.  A 
description  of  the  computer  program  contained  In  Appendix  A  for  use 
In  solving  maximal  flow  problems  for  large  networks  Is  Included. 

The  Generalized  Network 

Generally  maximal  flow  solutions  are  required  for  networks 
which  do  not  fit  the  restrictive  definition  of  a  basic  network.  For 
the  purpose  of  this  thesis  the  generalized  network  Is  defined  as  a 
modified  basic  network  with  a  set  of  arcs  A,  and  a  set  of  nodes  N, 
such  that: 

1.  Each  set  uf(1)  and  u>+(1)  can  contain  any  number  of 
directed  and  undirected  arcs. 

2.  Let  D  be  a  subset  of  A  such  that  If  and  only  If  the  arc 
a.j  j  Is  directed  It  Is  a  member  of  D. 

3.  Let  U  be  a  subset  of  A  such  that  If  and  only  If  the  arc 
a^  j  Is  undirected  It  Is  a  member  of  U. 

4.  K  Is  a  subset  of  N  such  that  If  the  node  1  Is  bounded 

It  Is  a  member  of  K.  To  the  bounded  node,  1,  we  associate 
two  real  numbers,  b^,  the  lower  bound,  and,  c^,  the  upper 

! 


bound  such  that  the  total  flow  through  the  node  1, 
denoted  y^ ,  Is  bounded: 

b1  ±  *1  i  C1  1  c  K 

5.  All  nodes,  j,  not  bounded  belong  to  the  subset  1,  where 
L  c  N. 

6.  The  network  may  or  may  not  contain  multiple  sources  or 
sinks,  but  shall  contain  at  least  one  source  and  one  sink. 

Procedure  for  Reduction  of  the  Generalized  Flow  Problem 

In  order  to  solve  for  the  maximal  flow  in  the  generalized  net¬ 
work  by  the  techniques  presented  In  Chapter  III,  the  network  must  be 
transformed  into  an  equivalent  basic  network. 

The  proposed  methods  for  reduction  are  based  upon  results  of 
the  general  theorems  presented  in  Chapter  II.  References  cited  give 
further  proof  of  some  of  the  methods  used. 

First,  consider  the  case  In  which  the  generalized  network,  G, 
contains  both  directed  and  undirected  arcs.  This  network  can  be 
simplified  by  replacing  each  undirected  arc  by  two  directed  arcs  as 
follows: 


Given  an  undirected 

1  arc 

ai  ,j  £ u 

Replace  a^  j  by  a^ 

j  and 

1  aj.i  ^  1 

Where 

Ci,j 

■  ci  .j 

Cj,1 

■  cfj 

Upon  reaching  a  solution  for  the  maximal  flow,  the  orientation 
by  the  flow  In  the  undirected  arcs  can  Immediately  be  determined  from 
the  values  of  y.  ^  and  y^  .  since  one  flow  will  be  zero. 
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Now  consider  a  network  in  which  two  distinct  nodes,  1,  j,  are 


connected  by  more  than  one  arc.  Successive  arcs,  flows,  and  bounds 
are  denoted  by  prime  symbols.  For  example  if  there  are  m  arcs  from 
node  1  to  node  j,  the  ktb  arc  would  be  denoted  a^  It  will  be 

assumed  that  each  arc  has  both  upper  and  lower  bounds. 

k  =  1 ,2,. .  .m 

■  »j  —  •  »j 

(k) 


b  (k>  <  y  (10  <  c  (k) 

D1.J  -y1.j  —  Ci » j 


The  set  of  arcs  a 


1 J 


k  *  (l,...m)  may  be  replaced  by  a  single  arc 

bj 


a  where*  c  s  0)  c  h  *  ?  h  (k) 

a1tj,  where,  ci,j  ci,j  b1,j  b1J 

Next  assume  a  pair  of  nodes  (1,j)  are  connected  by  two  oppositely 


directed  Ercs  a^  j  and  aj  ^  c  D. 

The  pair  of  arcs  may  be  replaced  by  a  single  directed  arc 
In  either  direction  as  follows:  If  c.  .  >  b,  . 

*  >J  J $  I 

Replace  a1  j  and  aj  1  by  a^  j  such  that:* 

C'l,j  =  C1,j  "  b1J 


b'u  =  b1,J  *  C1,j 

Otherwise  orient  the  arc  from  j  to  1. 

If  a  network  has  multiple  sources  and,  or  sinks,  merely  create 
a  single  source  and,  or  single  sink.  Then  construct  one  directed  arc 
of  Infinite  capacity  from  the  new  source  to  each  old  source  and  one 
directed  arc  of  Infinite  capacity  from  each  sink  to  the  new  sink  [8]. 

If  the  flow  through  the  nodes  Is  restricted  to  lie  between 
two  values,  the  character  of  the  maximum  flow  problem  Itself  remains 
unaltered.  Assume  we  are  given  a  network  In  which  the  flow  Is  subject 
to  both  arc  and  node  capacities,  where  all  arcs  are  directed  or  have 
been  converted  to  directed  arcs  by  the  previous  techniques.  This 
network  can  be  reduced  to  the  basic  network  maximum  flow  problem  by 
the  following  procedure.  Expand  the  network  such  that  for  each  bounded 
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re  created,  so  that 


node  (1)  c  K,  two  n»;v.  abounded  nodes  ( i ' ,  •«  *  * ; 
the  single  directed  arc,  a^ ,  j , ,  ,  joins  the  two  new  nodes  with  the 
upper  and  lower  bounds  associated  with  the  bounded  node  In  the  reduced 
network.  To  each  arc,  ak  y  create  ar.  arc  a'k  y  in  G'  and  to  each 
arc  a.  ,  In  G  create  an  arc  a.,,  .  in  G*.  The  bounds  on  the  other 

I  fj  '  *J 

arcs  remain  the  same. 


Example  Using  the  Simplified  Technique 

Given  the  ge.. trail  zed  network  showi  in  Figure  1,  we  are 
required  to  compute  the  maximal  flow  from  node  2  and  3  to  node  10. 


Or 


Maximize  F  =  T. 

y1  ,n  — 

The  flow  in  each  arc  is  bounded  as  follows 


(n)yi,n  "  y5,10  +  y7,10  +  y6 ,10 


3  -  y2.4  i8 

°iy2,6i4 

3  i  y3,4  i  7 

5  <  y3  5  <  10 

3  i  y4.5  i  8 

0  -y4,5  <-5 

(1) 

(2) 

0  iy4,6  -  3 

'  -  y4 ,6  -  2 

0  —  y4,6  —  2 

3  <  y5,4  i  8 

0  -  y5,6  —  8 

3  iy5,7i  7 

2  -y5,10±  5 

1  iy6,7i3 

3  iy6,10i  9 

(1) 

(2) 

1  —  y7,5  —  4 

1  1  y7,6  -  3 

2  -y7,6  -  5 

4  <y7  ioll0 

The  flow  In  all  nodes 

except  5  and  6  Is 

unrestricted. 

node  5,6  Is  restricted 

as  follows, 

4  <  y^  <  15 

4  <  yg  <  10 

This  flow  problem  can  be  simplified  by  using  the  techniques 
presented  In  the  last  section.  The  procedure  is  outlined  by  ordered 
steps.  Although  many  steps  could  be  combined  for  this  simple 
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FIGURE  1 

GENERALIZED  NETWORK  G  WITH  UPPER 
AND  LOWER  BOUNDS  ON  ARCS  AND  NODES 


15 

4 


10 

4 


FIGURE  2 


GENERALIZED  NETWORK  WITH  SINGLE 
ARCS  REPLACING  MULTIPLE  ARCS 


1  "■&****1t*a»mr* 


example,  It  Is  felt  that  +or  hand  and  computer  computation  this  Is  the 
most  efficient  procedure. 

1.  Add  the  single  sou^e  node,  denoted  as  1,  to  the  network. 

In  addition,  add  the  arcs  2  and  a-|  3  t0  the  network 
such  that: 

0  iy1>2  - +  “ 

0  1*1,3  -+  ” 

2.  Replace  each  undirected  arc  a^  j  with  upper  bound  j  >  0 
by  two  directed  arcs  a^  j  and  aj 

Therefore: 

a4>5  becomes  a^  and  a^J  0  <  <  5,  0  <  <  5 

131  f3l 

>4.6  becomes  >4,6  >"d  >6,4  *4,6  i2'  0  1*6,4  i  2 

>5.6  beCores  >5,6  >"d  >6.5  0  ±-*5.6  5  8-  0  i  *6,5  -  8 

3.  Combine  all  arcs  joining  each  pair  of  arcs  (1,j)  according 
to  the  following  procedure. 

a.  Combine  all  arcs  which  are  similarly  directed.  Into  a 
single  arc.  Thus  the  arcs  apg,  a42g,  a^  become  a4  g 


c  =  c{1) 

c4,6  c4,6 

+  C(2)  ♦ 

c4,6 

c(3) 

C4,6 

=  3 

+  2 

+  2  =  7 

b.  «  bP l 

4,6  4,6 

♦  til  ♦ 

4,6 

b(3) 

b4,6 

=  1 

+  0 

+  0  =  1 

Similarly  form  the  arcs  a4 

,5’ 

V 

\  and  a7,6  * 

C4,5  = 

13 

b4,5  =  3 

C5,4  = 

13 

b5,4  =  3 

C7 ,6  = 

8 

b7,6  =  8 

•**«  *!**■•  -a  >■  -  '■* . 


b.  Transform  each  pair  of  arcs  which  are  directed  In 
opposite  directions  between  each  pair  of  nodes 
Into  a  single  arc. 

Thus:  a^  g  and  ag  4  become  a^  5  3 

c4  5  *  c4  5  -  bg  4  =  13-3=10 

^4  5  *  ^  5  ‘  c5  4  *  3-13=-  10 

Similarly  the  pair  of  arcs  a^  g  and  ag  4*  a5  e  and 

a6  5*  a5  7  and  a7  5*  a7  6  and  a6  7  are  comb1ned-  The 
simplified  network  Is  shown  In  Figure  2. 

3.  Complete  the  transformation  to  the  basic  network  by  the 
following  procedure. 

a.  Expand  the  network  so  that  all  nodes  have  Infinite 
capacity  for  flow  as  outlined  earlier.  Thus  for 
nodes  5,  and  6  the  new  arcs  ag  g  and  ag  g  are 
created  and  the  equivalent  network  shown  In  Figure  3 
Is  established. 

b.  Add  the  return  arc  a-jg  ^  for  computational  purposes. 

4.  Transform  the  network  G  to  G'  by  using  the  change  of 


variables  equation 

yu  ay,u 

+  bi.j 

aU^A 

y10,l  s  y'l0,l 

yl  ,2  y'l,2 

yl  ,3  *  y  1 ,3 

y2,4  =y,2,4  +8 

y2,6  *  y2,6 

y3,4  =  y  3,4  +  3 

y3,5  =y'3,5  *  5 

y4,5  =  y4,5 

-  10 

y4,6  =  y  4,6  "  1 

y5,8  =  y  5,8  +  4 

y6 ,9  =  y6,9 

+  4 

y7,6  =  y*7,6 

y7,10  =  y7,10  +  4 

y8,6  =  y8,6 

-  8 

y8,7  =  y  8,7  "  1 

y8,10  '  y  8,10  +  2 

y9,10  =  y  9,10 

+  3 
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The  equivalent  problem  Is  obtained: 

Maximize  y'^Q  ^  with  the  node  flow  satisfying  the  balance 
equations  as  follows: 

Node 


1 

y’lO,l 

“  yl  ,3  y'l ,2 

-  0 

2 

y'l,2 

"  y'2,4  "  ^2,6  "  8 

*  0 

3 

y'l,3 

'  y'3,4  "  y3,5  ”  8 

-  0 

4 

y,3,4 

+  y‘2,4  ‘  y  4,5  '  y' 

4.6  +  22 

-  0 

5 

y  *3,5 

+  y,4,5  y5,8  "9 

-  0 

6 

y  2,6 

+  y,4,6  +y,7,6  +y' 

8.6  -^'6.9-  13 

=  0 

7 

y,8,7 

"  y ' 7 ,6  '  y,7,10  “  5 

a  0 

8 

y,5,8 

-y,8,6  "  y  8,7  -y' 

8,10  +  11 

*  0 

9 

y,6,9 

'  y  9,10  +  1 

»  0 

10 

y  *8,10 

+  y'7,10  +  y ' 9 ,1 0  "  y' 

10,1  +  9 

-  0 

Therefore  tne  values  of  the  remainder  are  Immediately 
seen  to  be  as  follows: 

d-j  *  0  dg  *  -  8  dg  *  -  8  d^  *  +  22 

dE  *  '  9  d6  "  '  13  d7  ’  *  5  d8  ’  +11  d9  ’  +1  d10  "  +  9 

Using  the  values  of  dj,  1  ■  1,,,,10,  as  the  capacities,  the 

new  source,  sink  and  arcs  are  added  to  obtain  the  equivalent 

network  G'  shown  In  Figure  4. 

5.  Then  using  the  Ford  Fulkerson  algorithm  to  obtain  a  feasible 
flow  In  G  the  following  values  are  obtained  for  the  flow  In 
G'. 
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y,1,2 

«  12 

y‘i.3-* 

y,2,4  “ 

0 

y2,6 

y  3,4 

=  0 

y3,5  *« 

y *4,5  * 

14 

y,4,6 

y,5,8 

■  5 

>'6,9  ■  6 

y  7,6  “ 

0 

y  7,10 

y,8,6 

=  6 

>'8,7  ■  7 

y  *  8 , 10  = 

3 

y9,10 

y ' 10,1 

*  20 

6.  Note  that  the  flow  in  all  sink  arcs  and  source  arcs  is  at 
arc  capacity,  therefore,  a  feasible  flow  exists.  Dropping 
nodes  0  and  11  plus  the  associated  arcs,  maximize  flow 
from  node  1  to  node  10  using  the  Ford  Fulkerson  algorithm. 
In  this  case  the  feasible  flow  computed  Is  the  maximum 
flow.  Then,  using  the  change  of  variables  equations  the 
resulting  maximal  flow  on  G  Is  determined.  Finally  G 
Is  simplified  by  eliminating  arcs  a5  g  and  afi  The 
maximal  flow  Is  20.  The  flow  In  each  arc  is  as  follows: 

>2,4  "8  >2,6'  4  >3,4"  3  >3,5  "5 

>4.5  "  7  >4,6  "  7  >5,4  ■  3  >5,7  "  7 

>5,10  '  5  >6,4  *  0  >6,5  "  2  >6,7  "  3 

>6.10  "  9  >7,5  '  1  >7,6  '  3  >7,10  "  6 

The  flow  In  node  5  Is  9  and  In  node  6  Is  9. 

Computer  Program  Description 

Since  It  Is  often  the  case  that  It  Is  necessary  to  solve  the 
maximal  flow  in  a  much  larger  and  more  complicated  network,  a  computer 
program  for  the  solution  of  the  generalized  maximal  flow  problem  as 
presented  In  this  thesis  has  been  written.  The  program  Is  written  In 
FORTRAN  IV  for  use  on  an  IBM  360  computer. 
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READ  IN 
NETWORK 


ALL  ARCS  INCIDENT 
TO  AND  INCIDENT 
FROM  THE  SAME  SET 
OF  NODES  ARE 
COMBINED 


NUMERICAL  VALUE 
OF  THE  TERMINAL 
NODE  IS  INCREASED 
TO  PERMIT  EXPAN¬ 
SION  OF  THE  NETWO! 


COMBINE  MULTIPLE 
ARCS  JOINING  TWO 
NODES  IN  OPPOSITE 
DIRECTIONS 


rRANSLATE  BOUNDED 
NODES  INTO 
BOUNDED  ARCS 


COMPUTE  CAPACITIE 
FOR  NEW  SOURCE  AN 
SINK  NODES  IN  G' 


CREATE  G' 


GENERAL  FLOWCHART 


CREATE  CHAIN 
FROM  SCUkCE  TO 
SINK 


LOAD  PATH 
ARRAY 


COMPOTE  MAXIMUM 
FLOW  INCREASE 


I  FLOW  IS  INCREASED! 
ALONG  PATH 


ELIMINATE  ALL 
|  ARCS  FROM  SOURCE 
AND  SINK  NODES 
OF  G' 


YES 


NO 


TRANSFORM  FLOW 
FROM  G'  TO  G 


(TRANSFORM  BOUNDED! 
ARCS  TO  BOUNDED 
NODES 


SEPARATE  APPRO¬ 
PRIATE  ARCS  AND| 
FLOWS  TO  TWO 
OPPOSITE  ARCS 
AND  FLOWS 


CHANGE  VALUE  OF 
TERMINAL  NODE  AND| 
COMPUTE  MAXIMUM 
FLOW  FOR  G' 


RINT  ARC  FLOWS 
AND  MAXIMUM 
FLOW  FOR  G 


GENERAL  FLOWCHART  CONTINUED 
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The  maximal  flow  In  generalized  networks  Is  computed  using 
the  methodology  outlined  In  this  chapter.  In  addition,  a  set  of 


routines  Is  added  to  the  program  which  will  permit  the  user  to 
transfer  several  networks  to  Input  data  d'rectly  without  first 
reducing  the  networks  to  basic  networks  by  hand  or  auxiliary  methods. 
However,  as  the  program  Is  written  each  network  must  have  a  single 
source  and  sink  and  all  undirected  arcs  are  to  be  represented  by 
two  oppositely  directed  arcs.  The  program  has  been  checked  out  and 
functions  properly. 

The  present  program  Is  written  for  a  50  node  network.  In 
order  to  solve  larger  networks  the  size  of  the  arrays  In  the  dimension 
statement  must  be  Increased.  The  general  flow  chart  of  the  entire 
program  is  given  on  pages  30  and  31.  For  the  Interested  reader  a 
detailed  flowchart  Is  presented  In  Appendix  A.  A  Fortran  listing 
of  the  program  Is  contained  In  Appendix  B  along  with  the  solution 
to  the  orevlous  example. 

Conclusions  and  Recommendations 

The  general  maximal  flow  problem  can  be  solved  more  efficiently 
by  use  of  the  method  developed  by  H.  Greenberg.  The  advantages  of  the 
technique  are  evident  when  solving  small  problems  by  hand.  The  tech¬ 
niques  presented  provide  a  systematic  approach  to  obtaining  a  solution. 
The  problem  of  computing  the  maximal  flow  In  a  mixed  network  *n  which 
the  flow  In  the  undirected  arcs  are  bounded  below  by  a  positive 
number  cannot  be  solved  by  this  method.  In  fact,  to  the  best  of  n\y 
knowledge  this  problem  is  unsolved  In  the  general  case  except  by 
complete  enumeration  techniques. 


32 


BIBLIOGRAPHY 


1.  Claude  Berne,  The  Theory  of  Graphs  and  Its  Applications,  John 
Wiley  and  Sons,  1958 

2.  C.  Berge  and  A.  Ghoul ls-Hourl ,  Programming,  Games  and  Transporta¬ 
tion  Networks,  John  Wiley  and  Sons,  1965. 

3.  R.  G.  Busaker,  T.  H.  Saaty,  Finite  Graphs  and  Networks,  McGraw- 
Hill,  1965. 

4.  G.  B.  Dantzlg,  Linear  Programming  and  Extensions,  Princeton 
University  Press,  1963. 

5.  G.  B.  Dantzlg,  D.  R.  Fulkerson,  On  Max  Flow  Min  Cut  Theorem 
of  Networks,  Rand  Research  Memorandum,  ft  M  -  T4T8  -  1,  1955. 

6.  G.  B.  Dantzlg,  D.  R.  Fulkerson,  Computation  of  Maximal  Flow 
In  Networks,  Rand  Research  Memorandum,  R  M  -  1489,  1955. 

7.  L.  R.  Ford  and  D.  R.  Fulkerson,  Flows  in  Networks,  Princeton 
University  Press,  1962. 

8.  D.  R.  Fulkerson,  A  Network  -  Flow  Feasibility  Theorem  and 
Combinatorial  Applications,  Canadian  Journal  of  Mathematics, 

Vol.  II  (1959) .  pp.  440-451 . 

9.  Gale,  Davlu,  A  Theorem  on  Flows  In  Networks,  Pacific  Journal  of 
Mathematics,  Vol.  7  (1957),  pp.  1073-1082. 

10.  Greenberg,  H.,  On  Generalized  Maximal  Flow  In  Networks,  unpublished. 

11.  Ford,  L.  R.  and  Fulkerson,  D.  R.,  Network  Flow  and  Systems  of 
Representatives,  Canadian  Journal  of  Mathematics,  Vol.  10  (1958), 
pp.  78-85. 


12.  ___,  A  Simple  Algorithm  for  Finding  Maximal  Network  Flows 

and  Applications  to  the  Hitchcock  Problem,  Canadian  Journal  of 
Mathematics,  Vol.  9  (1957),  pp.  210-218. 


33 


APPENDIX  A 


DETAILED  FLOW  CHARTS 

This  appendix  contains  the  detailed  flow  chart  of  the  computer 
program  in  addition  to  a  description  of  variables  ui:ed.  The  statements 
are  written  In  FORTRAN  IV  with  arrows  indicating  the  logical  flow 
through  the  system.  The  Input  and  output  statements  are  not  enclosed. 
Description  of  Variables 

Variable  Name  Description 

C(I,J)  Capacity  of  the  directed  arc  from  nr>de  i  to  node  J. 

B(I,J)  Lower  bound  of  the  dii-ected  arc  from  I  to  node  J. 

Y(I,J)  Initially  used  to  read  In  data. 

Y(I,1)  The  number  of  che  node  the  arc  a.  .  Is  Incident 

from.  1 


Y( I ,2)  The  number  of  the  node  the  arc  a^  j  Is  Incident  to. 

Y(I,3)  Upper  bound  of  arc  a^  j. 

Y(I,4)  Lower  bound  of  arc  a^  y 

Y(I,0)  The  value  of  the  flow  In  a  directed  arc  from  I 
to  J  after  Input  data  Is  transmitted  to  arrays 
C(I,J)  and  B(I,J). 

XL ( I )  Dumny  array  used  to  compute  values  of  the 

capacities  of  the  source  and  sink  arcs  for  G' 
and  then  the  label  of  the  nodes  In  the  label 
routine. 


NL(I)  Dimry  array  used  in  the  label  routine  as  a  check 

to  determine  If  all  nodes  have  been  checked  and 
then  to  store  the  path  used  to  Increase  flow. 

XNC(  1,1)  Number  of  the  node  that  is  bounded. 

XNC(I,2)  Upper  bound  on  a  flow  through  a  node. 

XNC(I,3)  Lower  bound  on  flow  through  a  node. 


XNC(I,4)  Value  of  flow  through  a  node. 

NN  Number  of  networks. 


N  Number  of  nodes  In  the  original  problem. 

M  Value  of  the  designated  number  of  the  slnh  node. 

KX  Value  of  the  designated  number  of  the  source  node. 

NA  Number  of  arcs  In  th*  original  problem. 

NC  Number  of  bounded  nodes. 

Input  Card  Format 

All  entries  on  the  card  format  In  this  section  are  right, 
justified  and  are  fixed  point  entries  if  the  variable  r^me  begins 
with  I,  J,  K,  L,  M,  or  N  and  are  floating  point  entries  for  all  other 
variable  names. 

Card  #1 

Column  Name 

1-4  NN 

Card  #2 

1-4  N 

4-8  NA 

8-12  NC 

Card  #3  through  Card  NA  +  1 

1  -  10  Y{ I ,1 ) 

11  -  20  Y(I ,2) 

21  -  30  Y(I ,3) 

31  -  40  Y(I ,4) 
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Card  #NA  +  3  through  Card  NC  +  3 


Column  Name 

1  -  10  XMC(1 ,1 ) 

11  -  20  XNC(I,2) 

21  -  30  XNC ( I ,3) 

Detailed  Flowchart 

The  detailed  flowchart  of  the  program  Is  given  In  the  next 
11  pages. 


DETAILED  FLOWCHART 
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C(K,J)  =  -  C(K,J) 


C(L,J)  =  C(K,J) 
B(L.J)  «  B(K,J) 


DETAILED  FLOWCHART 


DETAILED  FLOWCHART 
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C(M-1 ,2)  -  .9999 
MX  -  M-l 


DETAILED  FLOWCHART 
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detailed  flowchart 
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PRINT 
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DETAILED  aOWCHART 


APPENDIX  B 


FORTRAN  -  IV  LISTING 

This  appendix  contains  the  complete  lasting  of  the  computer 
program.  The  problem  niustrated  In  Chapter  IV  was  used  to  display 
the  output  format.  Comment  cards  are  Included  for  clarification  and 
Identification  of  the  main  sections  of  the  program. 
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DIMENSION  C( 50, 501  ,8(50,501 ,Y( 50, 50), XL (50 1 ,NL( 501, 
lXNCf 50,4) 

READ ( 5.801 )  NN 
801  FORMAT < I  4) 

DO  800  I U=1 » NN 
READ (5,1)  N,NA  ,NC 

1  FORMAT ( 314 ' 

00  21*1. NA 

RE  AD ( 5 ,3 ) (  Y(I.J),J*1,4) 

3  F0RMAT(4F10.4  ) 

2  CONTINUE 

00  8  I=1,NC 

READ ( 5.4 )(XNC(  I,  J) ,J*1,3) 

4  FORMAT  *  3F10.4) 

8  CONTINUE 

KX*1 

M=  N  +  NC*  2 
RET*0 

DO  5  1=1,50 
NL  (  I  )  *0 
XL (I)  =0 
DO  5  J=1 ,50 
B(  I  ,  J)=0 
C( I, J)=0 

5  CONTINUE 

C  THE  CAPICITY  MATRIX  IS  LOADED  AND  ALL  ARCS  INCIDENT  TO 
C  AND  INCIDENT  FROM  THE  SAME  NODFS  ARE  COMBINED 
00  6  1=1. NA 
L=  Y ( 1,11  ♦  1. 

K=  Y(1, 2)  +  1. 

C(L,K)=  Yd, 31+  C(  L,K  ) 

B(L,K)=  Yd,  4)*  B(  L  ,  K  ) 
b  CONTINUE 

'•  _  NUMBER  I  CAL  VALUE  OF  THE  TERMINAL  NODE  IS  CHANGED  TO 
hLLOW  FOR  EXPANSION  OF  THE  NETWORK 
L=N+1 
K-N  +  NCH 
DO  11  1  =  1. M 
C(  I,K)=C(f,L) 

8(I,KJ=B(I,U 
C( I,L)«0. 

11  B(  I  ,  L )  =0. 

C  COMBINE  MULTIPLE  APCS  JOINING  TWO  NODES  IN  OPPOSITE 
C  DIRECTIONS. 

DO  15  1*1, M 
CO  15  J*1»M 

VF(C(I,J).GT.  O..AND.  C(J,I).GT.  0.)  GO  TO  16 
GO  TO  15 

16  It~(C(I,J).LE.B(  J.I  )  I  GO  TO  17 
C! I, J)*  C ( I , J ) -B ( J , I ) 

Bd,  J)=  B(  I  ,  J)  -C  (  J  ,  I  ) 

c(j,n*-c(j,:‘ 

GO  TO  15 

17  C( J, I)=C(J,I)-  8(1, J) 

B(  J,  I)=B(J,D-  C(I  ,J) 
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1 


C(I, J)=-C< I 
CONTINUE 


J) 


i 


156 


15  _  .  . 

TRANSLATE  BOUNOEO  NODES  INTO  BOUNDED  ARCS 
IF(NC) 155* 155* 156 
L*N*  1 

DO  34  1*1, NC 
K*  XNC (  I  .  1 )  + 1 . 

DO  32  J*l.  M 
IF(C(K, J) )  32.32,33 
C<L,J)=  C(K,J) 

8( L, J) *  B<  K ,  J) 

IFICIK.J)  .GT.O.I  C C K » J) *-C f  K,  J I 
CONTINUE 


33 


32 


C(K,L»*  XNC (  1,21 
B(K,Ll*XNC( 1 *31 


L-L+i 

34  CONTINUE 
155  CONTINUE 

C  THE  FOLLOWING  OPERATION  COMPUTES  THE 
C  ARCS  FROM  AND  TO  THE  SOURCE  AND  SINK 
DO  9  1*1, M 
SI*0 
S0*0 

DO  10  J=1,M 
Y( I, J1*0 

IF (C ( I , J) . GT .  Q.)  CC I ,  J )  *  C ! I , J) -Bf I , J ) 
IF(CU,J). GT.O.I  SO»SO*BU.J) 

IF f C ( J.II.GT.O)  SI*SI+B(  I-  I ) 

CONTINUE 


CAPACITIES  OF 
NODES  IN  G* 


THE 


10 


13 


4 

2 


XLII)-SI-SO 
9  CONTINUE 

G'  IS  COMPLETED  BY  ADDING  NEW  SINK  AND  SOUPXE  ARCS 
CIM-l, 21-9999 
MX*  M-t 
00  12  1*2,  MX 
IF  f  XL ( I ) I  13.12,14 
CCI.Ml*  -XL  111 
GO  TO  12 

cu.n*  xliii 

CONTINUE 
00  31  1*1,50 
XL  ( I  1*0 
XL  (  KX  )  — K X 

LABEL  ROUTINUE  AS  PRESENTED 
USEO  TO  DETERMINE  THE  FLOW  CAN 
00^210  L*KX  ,  M 

IF  I  XL (  1 1  193,93,96 
NL  ( I  )*0 
GO  TO  22 
LN*L 

DO  94  K*l.LN 
IF ( NLI K  J 194,94,95 
I*K 

GO  TO  96 


3l 


2 

THE 

IS 


BY 

BE 


FORD  AND  FULKERSON 
INCREASED 


96 


93 


95 


50 


:A  ,  \k' "  W  'T 


-Try 


%  v*v  ■  r  at. 


c 

c 


94  CONTINUE 
GO  TO  210 

22  00  21  J=  KX f  M 
IF5XLI J) >23,23*21 

23  IFICII.J))  24,24,25 

25  IF  l  Y(I,J).LT.  C(I,JM  XL  C  J  >  =  I 
NJ.C  JI«XL  ( J) 

GO  TO  21 

24  IF (C I J , 1)121,21,26 
*  IF  j  Y( J, I ).GT.  0.  IXL( J)>I 

NL  i  J )  =) . 


26 


30 


'XL  ( J ) 

21  CONTINUE 
GO  TO  93 
210  CONTINUE 

IF  (  XL  I M) 127,27,28 
A  PATH  FROM  SOURCE  TO 
FROM  THE  LABELLED  NODES 
28  1  =  1 

NL ( I ) =M 
K=NL ( I ) 

1  =  1  +  1 

NL ( I )=XL I K) 

K=NL ( I ) 

IF (  K  .EQ.  KXIGO  TO  40 
GO  TO  30 
4C  X=9999. 9 

THE  MAXIMUM  AMOUNT  THAT  THE  FLOW  CAN  RE 
ALONG  THE  PATH  IS  COMPUTED 
45  DO  41  1=1, M 
K=NL ( I ) 

L=NL ( 1+ 1 ) 

IF (L ) 50, 50,  42 

,L))43,43,44 


SINK  IS  DETERMINED  DIRECTLY 


INCREASED 


42  I F  ( 


43  IF(fc(L,K»-  Y(L»KI).LT.  X»  X=  C  (  L , K ) ~Y ( L , K  ) 
GO  TO  41 

Y(K,L> 


XI  X' 


44  IF (  YIK.L1.LT. 

41  CONTINUE 
50  CONTINUE 
THE  FLOW  IS  INCREASED 
DO  51  1=1, M 
K=NL 1 1) 

L=NL( I+l 1 
IF (L ) 60, 60,  52 

52  IF(C(K,UI  53,53,54 

53  YIL.K1*  Y( L , K ) ♦  X 
GO  TO  51 


54  YI K,L ) =  Y( K , L  J  -X 
51  CONTINUE 

60  CONTINUE 

DO  61  1=  KX , M 
NL ( I )=0 

61  XL ( I )=0 
GO  TO  20 

C  i  HE  SINK  CANNOT  BE  LABELLEO  AND  A  CHECK  IS  MADE 
C  TO  GERMING  IF  THE  FLOW  IS  FEASIBLE 


27  IF  C  M.EQ. (N*NC+1) I  GO  TO  80 
MJ-M-1 

DO  70  1^  2,  MJ 

IFdCCl.II-  Ydtin  .GT.  .0001  .OR.  (C(l»M|- 
lYd.MIl.GT.  .0001)  GO  TO  90 
70  CONTINUE 
KX»2 

1 

CIM.KXI-0 
GO  TO  20 

C  THE  ORIGINAL  NETWORK  IS  RECONSTRUCTED  AND  FLOW  IS 
C  TRANSFERED  FROM  G*  TO  G 

80  CONTINUE 

DO  81  I *KX  ,  M 
DO  81  J«KX,M 

IFCCCI  , JI.LT.O.)  GG  TO  81 
CCI,J)»  CU.JI  +B(  I  •  J ) 

YU,  J)»yTi,J)*B(I,J) 

81  CONTINUE 

IF(NC) 157,157,158 
158  L»N+1 

DO  100  I  -1 » NC 
K»  X  NC  (  I  ,  1 1  + 1 . 

XNC< 1,4)  *  Y(K,L) 

C( K, L)«0. 

YCK,L)«0 


M  Li 

DO  101  J-KX.M 

Tf(cTl,j).le.i 

C(K,J)=C(L,J) 
Yf  K, J)*Y(L , J) 
Y(L,J)»0 
CCL, J)«0 


0)  GO  TO  101 


ioo  Continue 


157  CONTINUE 

DO  102  I«  KX,M 
DO  102  J-  KX.M 
IFCCC J, I ))  103,102, 102 

103  Cl  J,  II— Cl  J,I 

C( I, J)*  CII.J  *B(  J, I ) 

B(  ,J)«  BCI.J  ♦CCJtU 

IFC ( B( I.J)-B(J,I) ) .GT.  Y ( I , J 1 1  GO  TO  104 
YC J, I)»6(J,I) 

YU,  J)»Y(I,J)  +  B(J,I) 

GO  TO  102 

104  Y(j,n*B(I,j|-Y(I,J| 

Y( I, J)*B( I , J) 

102  CONTINUE 
L»N+1 
K-N+NCU 
DO  105  I *KX , M 
C> I,L)«CCI,K) 

B( I ,L)*B( I ,K) 

YC I,L)*YCI,K) 

105  C( I , K) *0 • 


r 


«n»y 


A* 


WRITE(6.86) 

8fc  FORMAT ( lHl t  25X,  •  ARC • , 9X , • LOWER  BOUND' , 12X ,» FL DM • , 
U4X. 'UPPER  BOUND',//) 

DO  82  1=2, M 
RET=RET+Yt I ,M| 

DO  82  J=2,M 
I F  C  C  C I ,J) >82,82,85 
85  L*I- 1 
K=  J- 1 


WRITE (6. 83)  L,K,B!I,J),Y(I,J),C(l,J) 

8?  F0RMAT(24X,I2,1X,I2,  5X ,F1 0. 3 , l IX, F10. 3, 10X , F 10. 3, / ) 
82  CONTINUE 
WRITE(6,87) 

87  FORMAT!  1HO  ,24X»  'NODE'  ,9X,  'LOWER  BOUND'  ,  12X,  •  FLOW  . 
114X, 'UPPER  BOUND',//) 

IF(NC) 161 ,161,160 

160  DO  89  1=1, NC 
K=XNC (1,1) 

HR  ITE ( 6 , 88 )  K  ,XNC(I,3),XNC(I,4), XNC  (1,2) 

88  FORMAT(25X,I2,6X»F10.3»llX,FlO.3»lOX,F10.3»//) 

85  CONTINUE 

161  CONTINUE 

MR ITE ( 6, 84 )  RET 

84  FORMAT (  25X,'MAX  FLOW  =',F10.4) 

GO  TO  91 

90  WRITE  (6,92) 

92  FORMATI25X, 'NETWORK  HAS  NO  FEASIBLE  FLOW') 

91  CONTINUE 
DO  1J7  1=1,50 


00 


Y(  I.J)=0 
117  CONTINUE 
800  CONTINUE 
END 


7  J=1 , 50 
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